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ABSTRACT 


An investigation related to the use of linear indicial 
functions in the time and frequency domains, enabling 
one to derive the proper aerodynamic loads as to study 
the subcritical response and flutter of swept lifting 
surfaces, respectively, of the open/closed loop 
aeroelastic system is presented. The expressions of the 
lift and aerodynamic moment in the frequency domain 
are given in terms of the Theodorsen’s function, while, 
in the time domain, these are obtained directly with the 
help of the Wagner’s function. Closed form solutions of 
aerodynamic derivatives are obtained, graphical 
representations are supplied and conclusions and 
prospects for further developments are outlined. 

NOMENCLATURE 

a n Dimensionless elastic axis position measured 
from the midchord, positive aft 
c n Chord length of wing, normal to the elastic axis, 
2b, 

C L(Xn Lift-curve slope 

fh i fa Plunging and pitching deflection functions 
h , fy, Plunging displacement and its amplitude, 
respectively 

I a Mass moment of inertia per unit length of wing 

/ Wing semi-span measured along the midchord 
line 

l h Dimensionless aerodynamic lift, L^/mU* 

Aerospace Engineer, Ph.D. Student. 
f Professor of Aeronautical and Mechanical 
Engineering, Department of Engineering Science and 
Mechanics. 

* Senior Research Scientist, Senior Aerospace 
Engineer, Aeroelasticity Branch, Structures Division, 
Senior Member AIAA. 

Copyright © 2000 The American Institute of 
Aeronautics and Astronautics Inc. All right reserved. 


L b ,l b Overpressure signature of the N-wave shock 
pulse and its dimensionless counterpart, 

L b bjmU 1 2 , respectively 
m Airfoil mass per unit length 

m a Dimensionless aerodynamic moment, M a b 2 / I a U 2 
N Load Factor, h'/g 

P m ,p m Peak reflected pressure in excess and its 

dimensionless value P m b/mUl , respectively 
r Shock pulse length factor 

r a Dimensionless radius of gyration, (l a /mb 2 } ~ 

5 ,£ Laplace transform variable and operator 
S a ,x a Static unbalance about the elastic axis and its 
dimensionless counterpart, S a /mb 
t , t {) Time variables 

U oot U n Freestream speed and its component normal to 
the elastic axis 

V Dimensionless free-stream speed, U Oa jbo) a 

x Coordinate parallel to freestream direction 

x Chordwise coordinate normal to the elastic axis 
y Coordinate perpendicular in the freestreem direction 
y Spanwise coordinate along the elastic axis 
w Downwash velocity 

z Transverse normal coordinate to the 

midplane of the wing 
Z Vertical displacement in z direction 

a , a {] Twist angle about the pitch axis and its 
amplitude, respectively 

£/i’?a Structural damping ratio in plunging, c h j2mo) h 
and in pitching, c a /2I a 0) a 

77 Dimensionless coordinate along the wing span, y / 1 
X Spanwise rate of change of twist 

A Swept angle (positive for swept back) 

/i Airfoil/air mass ratio, mjizpb 2 

^ Dimensionless plunge coordinate, hlb 
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p Air density 

o Spanwise rate of change of bending 

T Dimensionless time, U n t/b n 

r p Dimensionless positive phase duration of the 
pulse, measured from the time of the arrival 
co,k n Circular and reduced frequencies, cob n /U ft 

Uncoupled frequency in plunging, {K h /m) U2 

and pitching, {K a /I a )' n 
CO Plunging-pitching frequency ratio, co h /co a 

Subscript 

( • ) c Circulatory terms of lift and aerodynamic 
moment 

( * )„ c Non-circulatory terms of lift and aerodynamic 
moment 

( • )„ Quantity normal to the elastic axis 
A ( • ) Quantity associated with the swept wing 

Superscript 

( • ) Variables in Laplace transformed space 
(•),(•) Derivatives with respect to the timer , and 
the dimensionless timer , respectively 

INTRODUCTION 


In this paper, the linear indicial functions in the time 
and frequency domains are used to determine the 
corresponding unsteady aerodynamic derivatives for 
swept lifting surfaces. Such a treatment of the problem 
enables one to approach either the open/closed loop 
aeroelastic response in the subcritical flight speed 
regime to arbitrary time-dependent external excitations 
(such as e.g. gusts, airblasts due to explosions or sonic- 
booms), or the flutter instability of actively 
controlled/uncontrolled swept wings. Moreover, this 
study is intended to extend its scope as to include, by an 
analogous procedure, the case of the various flight 
speed regimes, i.e. in addition to the incompressible 
one, also the compressible, transonic (within the 
linearized concept of indicial functions), and 
supersonic. 

The representations of the aeroelastic governing 
equations in the time-domain is useful toward the 
approach of both the subcritical aeroelastic response of 
the actively controlled/uncontrolled swept aircraft 
wings. In other worlds, in this framework, the 


open/closed loops dynamic response of aeroelastic 
systems can be analyzed. 

The unsteady aerodynamic lift and moment in 
incompressible flight speed regime are expressed for 
the swept aircraft wing in the time and the frequency 
domains by using the Wagner and Theodorsen 
functions, respectively. For the response of dynamic 
systems it is only necessary to express the lift and 
moment via the indicial Wagner’s function. For the 
approach of the flutter problem, the Theodorsen's 
function helps the conversion of the expressions of both 
the aerodynamic loads and the unsteady aerodynamic 
derivatives in the frequency domain. Herein the case of 
a 2-D lifting surface, including the plunging and 
pitching degrees of freedoms is considered. 

PRELIMINARIES 


In the next developments an extensive use of variables 
in both the time and frequency domains will occur. As 
shown in Edwards, Ashley & Breakwell 2 3 , for zero 
initial conditions, the aeroelastic equations can be 
converted from the time to the frequency domain via a 
Laplace transform. This results in the possibility of 
using the correspondence s — > ik n , where s and k n are 
the Laplace variable and the reduced frequency, 
respectively. The Laplace transform operator £ is 
defined as: 


(i) 

In this sense, the Wagner's function 0( r) is connected 
with Theodorsen’s, function C(^) via Laplace 
transform as: 


C(fc„) _ *•(*„)+ «G(0 

<k„ ik„ 

and vice-versa: 


£>(r)r'*" r t/T = <d(;/c„ ), (2) 


<Kt )=£‘{c{k n )lik n }, Re(ik n )>0. (3) 


and having in view the correspondence 5 ik n we can 
also write: 


®{ik„) = £>(r> _sr <ir = 0(s), (4) 

ik n —>s 

Using this relationship, it is possible to obtain the full 
expression of unsteady aerodynamic coefficients in 
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terms of the Theodorsen's function C(k n ) and its 
circulatory components F(k n ) and G(k n ). 

It is interesting to note that the reduced frequency 
parameter k n for swept and for straight wings coincide: 


where h = h(yj), a=a(y\t) are the displacements in 
plunging and pitching, respectively, and the origin of 
the x axis coincides with the elastic center, the 
downwash velocity vv normal to the lifting surface 
becomes: 


£ _ _ col? cos A _ cob 

U n U cos A U 
icot = ik n T . 

This implies that the indicial Wagner's function <p(r) 
remains invariant to any change of the sweep angle. 

ANALYTICAL DEVELOPMENTS 


( 5 ) 

( 6 ) 


w(jc, y,t)= vv(jc, >’,?)= dZ/dt + U „ cfLjdx . ( 11 ) 

In conjunction with the definition of Z, (Eq. 10), Eq. 
(11) becomes: 

w(x, y,r)= h + xa + U^a cos A 

+ U „{dh!dy + xdajdy) sin A . (12) 


For swept wings, the total lift per unit span, can be 
expressed in the form: 

Ak i G’’ * ) 

+ ^ndCy,ty^ nc - i {y,t)- ( 7 ) 

where the indices c and tic identify the various 
contributions associated with the circulatory and 
noncirculatory terms respectively. 

Using similar notations, the total moment per unit span 
about the elastic axis is: 

fi M a (y,t}= A M c {y,ty^M nci {y,t) 

+*M ncl {y,t)¥ K M„'i(y,t)± A M a {y,t), ( 8 ) 

A M a (y,t ) being associated with the apparent moment 

of inertia. Herein the lift is positive in the negative z 
direction (considered positive downward), while the 
moment is positive nose up. For the sake of 
convenience, herein the plunging coordinate is positive 
when is downward (see Fig. 2). 

Using the expression of the lift (Eq. (7)), the equation 
for the moment, Eq. (8), can be cast as: 


where the superposed dots denote the derivatives with 
respect to time t. The quantity in Eq. (12) undescored 
by a solid line is usually discarded in the specialized 
literature (see Bisplinghff, Ashley & Halfman 2 ), as 
being related with the wing camber effect (see also 
Flax 4 ). However, herein this quantity will be included. 
The in-plane coordinate x normal to the elastic axis 
(see Fig. 1) can be expressed as: 


x = b„(\/2-a n ). 03) 

Consequently, using the dimensionless timer(= U n tjb n ) 
Eq. (12) becomes: 


Xx,y,T)=U n 


h' 


dh 


+ a + — tan A 
b„ dy 


da 




+ a' + ^— tanA j | (14) 


where (•)= d{)/dr . 


In the following sections, the unsteady aerodynamic 
loads in incompressible flow can be obtained in time 
and, with the use of the Laplace transform space, in the 
frequency domain. 


A M a (y,t)=-{\/2 + a„)b„ A L c (y,t)~ a n b„ A L nci (y,r) 

+ 0/2 - a n )b n A L„ c 2 (y, t >- A M „ c3 (y, t) 
+ A W„(^0. (9) 


Unsteady Aerodynamic Loads in Incompressible 
Flow 

Time Domain 


Basic Considerations 

Expressing the vertical displacement Z of a point on the 
center line of the wing as, (Fig. 2): 

z(x, y, t)= h + xa , (10) 


The circulatory components of the lift and moment 
expressed in term of Wagner’s indicial function <p(r) 
(called also heredity function) obtained in the time 
domain are: 
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X (y>T ) = ~ C La, b n P U n L 0 ( T ~ T <> ) — ^ ~ <ft<> * 




lL A# e (y,x) = 


1 

— + fl„ 
2 " 


fLaMpUn 




dr n 


(15) 


(16) 


As concern the aerodynamic non-circulatory 
components, using the dimensionless time r these are 
expressed as: 


a ^ (3^)= -1 C tan pU 2 „ [h' - a„b n a’l (17 a) 

,L nc2 (y,t)=~C Lan pU 2 n by, 07 b) 


i L „c 3 (y>'r)=--C Lan pU 2 n b; tan A 


X+l)f + 5 r A 


+s r 

+ 5 r 


<9tr 

tan A 

<9y 

t?A 

— tan A . 

# J 


+-c LanP uy n 


tan A 


( ^ +l) f 

(17c) 


The moments induced by the non-circulatory 
components of lift are expressed using the 
dimensionless time as: 


i M ncl (y.T )=jC Lan pU ; [h' - a„b n a'\i„b „ , (18 a) 


M nc Ay,T)=-'-C Lan pUy n (l 


2 a " 


P . 


(18 b) 


a M „ c3 (y,r)=-ic Lan pUy^X tan A 


+ tan A| 


c <9<r 

+ <5 r — tan A 
dy 


1 


(<5, + l)^ + S r A 
b n 


c La pvX\-+< 


8 


tan A 


(5 r +l)A + 5 ^tenA 

3y 


X a (y^)=~pc La yuy. 


, (18 c) 
(18 d) 


as in the following ones, the terms affected by the tracer 
S r identify those generated by the last term in the 
expression of the downwash velocity (Eq. (12)), (term 
underscored by a solid line). In general these terms are 
discarded being considered negligibly small (see 
Bisplinghoff et al.\ and Flax 4 ), in what case S r = 0 , 

otherwise 8 r = 1 . Replacement of Eqs. (15) and (17) in 
Eq. (7) and of Eqs. (16) and (18) into Eq. (8), results in 
the unsteady lift and aerodynamic moment expressed in 
the time domain. Concerning the circulatory parts of the 
lift and aerodynamic moment, Eqs. (15) and (16), these 
can be explicitly determined by transforming these 
expressions in the Laplace domain using the 
relationship between the Laplace transform of the 
Wagner and Theodorsen functions, namely 
C(- is)/s = Jf(^(r)) = O(s), and afterwards inverting 
back the obtained expressions in the temporal space. 
Alternatively, in order to ease the computations, the 
available approximate expressions for <j>(j) (see 

Edwards et al. 3 ) can be used in the Laplace transform 
process. 

The expressions of lift and aerodynamic moment in the 
time domain, A L h (y y r) and A A7 a (y,r) can be used to 
determine the subcritical aeroelastic response of swept 
wings. However, when the aeroelastic response of 
lifting surface to time-dependent external pulses, is 
needed, the unsteady aerodynamic loads in the time 
domain, A L h and A A/ a , have to be supplemented by 

the ones corresponding to above mentioned pulses. 

This will be considered in the next developments and a 
simplified illustration of the capability of this method 
will be given in this work. 

Laplace Transformed Space 

Several preliminaries related to Laplace transform 
applied to aeroelastic quantities will be given next. 

=lA* L Ay’A*M a {y,T)y sT dT . (i9) 

We will express the Wagner’s function and the 
plunging and pitching degrees of freedom in the 
Laplace transformed space as: 

<&{s) h(t) h(s) <*(y) 


where the spanwise rates of change of bending and Considering zero initial condition, Laplace Transformed 

twist, a and A, respectively, are expressed as counterparts of Eqs. (15) and (16) are: 

O-dhjdy and k-dajdy . In these equations as well 
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r 


2 h 

s — + sa + so tanA + 
b 


(23 d) 


L,{y^-c L hpu: 


|s 2 d + fc„sAtanA)]o(j) ,(20) 


,M c (y,s)= 


(\ 


— + a n 
2 


f La b nP U 2 n 


(1 

— ~ a ~ Itf 

l 2 


2 h 

s — + sa + s<7tanA 


2 a + b n sX tan a) p(.s) . (21) 


Herein, and in the next expressions, the plunging and 
pitching motions in the Laplace domain are expressed 

as h = h(y, 5 ) and a = a(y, s) . 


The Laplace transformed counterpart of Eqs. (17) can 
be expressed as: 


l Kc\ (y> S ) = “ C La n P U n s2 L /2 - a n b n& l ( 22 a ) 


i L nc 2 (y^)=-^c Lan pu'-b n sd , 
,L nci {y,s)=-U Lan pUy n i^A 


(22 b) 


(l + <5,)f * + S r A 
b n 


~ do 

+ S r — - tan A 
dy 

+ 8 r — tan A 

dy 


+ -CLa n P U n b l tan A | 


0 + *,)f* 

b„ 


(22 c) 


The Laplace transformed counterpart of Eqs. (18) can 
be written as: 


A M „ fl (y. *) = I C ta> pU 2 s 2 [5 - « a4 A , (23 a) 


,M nc2 (y,s)=--C Lan pUy n \--a n U, (23b) 


,M nc3 {y,s)=-±C Lan pUX 1 A tan A 


+ ^Cu,'PUya„ tan Aj 


(<5 f +l)— i + 5,A + 
b„ 


e <?<7 

+ 5 r tan A 

dy 


-\c LarP uy 




8 n 


(5 + 1 )— i + 8 tan A 
b„ ' 3y 


tan A 


(23 c) 


,M a (y,s)=~pC La y n u 2 n s 2 d. 

The equations of lift and aerodynamic moment help us 
to perform the conversion in the frequency domain. 
This will be done in the next section. 

Frequency Domain 

Upon replacing s — > ik n in Eqs. (22) and (23); using the 
relationship between Laplace transform of Wagner and 
Theodorsen’s functions (Eq. (2)); representing the time 


dependence of displacement quantities as: 

«(>’.d= fa (>’)«L *„)=/« (>>(/‘” r > 

(24 a) 

h(y\r)= f h (j)ft (t, k„)=f h {y]h 0 e ,k,,r , 

(24 b) 

and expressing: 

,L h (y,k n ,t)= A L h (y,k„y k -\ 

(25 a) 

a (y,k n ,r)=iM a (y,k„ )e'*” r , 

(25 b) 


the equations for the unsteady lift and moment 
amplitudes can be expressed in the frequency domain. 
These expressions that coincide with the ones obtained 
differently in Barmby et al.\ can be used in the flutter 
analysis of swept aircraft wings. 

In this analysis f a (y) and f h {y) are chosen to be the 
decoupled eigenmodes in plunging and twisting of the 
counterpart structure, and are determined as to fulfill 
identically the boundary conditions. Using the spanwise 
dimensionless coordinate rj = y/l, these are expressed 


as: 


fh(y) = Fh( 7 l) = C\ ((cos 


Pfl - cosh /?,77) 


sinh f3 ] + sin /3, 
cosh /Jj + cos /?, 


-h sinh /?j77 - sin fitf) 
fa(y)= F a( 1 l)= C 2 sin PlV . 


(26) 

(27) 


where for the first bending and torsion we have 
/?, = 0.5969;r and p 2 =/r/2 . 

The constants C, and C 2 are chosen as to normalize 
f h (y) and f a {y), and so to get the unitary maximum 
deflection at the wing tip. The uncoupled first bending 
and torsion mode shapes needed for the evaluations of 
the aerodynamic lift and moment are shown in Fig. 3. 
The following expressions will be useful in the next 
developments: 


A' = 


d 2 a 


3ydx 3y 


a n ik.e ik * r , 


(28) 
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W" 


ff /_ _ ffuGO 
<fcr d 2 f h (y), _«., 

~W‘~ h " e ■ 


(29) where, for the sake of convenience, the unsteady 
aerodynamic derivative are written as: 


(30) 

(31) 


//, = = //,;//, = (//, -// 5 )ff 4 = (» 4 

A, = A, ; A 2 = A, ; A, = (A, - A, > A 4 = (A 4 - A 6 ). (38) 


The unsteady aerodynamic derivatives in the frequency 
domain for swept wing will be obtained from the 
following equations, expressed in terms of Wagner’s 
function <&(ik n ). Comparing Eqs. (36) with that 
At this point, a careful inspection of Eqs. (7) - (8), expressing the lift in the frequency domain, yields: 


Unsteady Aerodynamic Derivatives in the 
Frequency Domain 


suggests the following representations for the lift and 
aerodynamic moment: 


ti 


. \ ^ r 1 C La„ , 

L \ =T C laJh -T— rA tanA 


dy 


L h (y,k n ,x)=-pU 2 „2b„ k„H ] — + k n H 2 a' + k 2 H 3 a + <5 A 


ds’ 1 


tanA 


+ <u(/* n f k n f h - ib n % tan A 


+ k 2 H, — + H,a'+H — 
b. b, t 


(32) 


dy 

(39 a) 


L , = — 


,M a (y,k n ,T)=±pU;2b; 


( 


h! 


-p-H'*.) if a !*/./<* an A 


k „\-r + k - n A. 2 a' + k;.%a 




, 2 . h . , . h* 

+ k n A A —+ A 5 a + A a — 

b„ j 


(33) 


(5 r +l)-^-A-„ +£,£„ — ^ tan A 
-A <?y 2 

1 Q 


1 b 

+^ C La n a n 77 tan A 

1 k n 

c 1 n b n df a i , -j.a 

-o — C /a — — —tanA — C lfl a„f n — f n . 

r i La " k 2 dy ° ^ J a 


2 JL 


(39 b) 


Herein the dimensionless unsteady aerodynamic From a similar analysis, comparing the Eqs. (39)'with 

coefficients H t , A i have been introduced, and k„ has that expressing the aerodynamic moment in frequency 


been included as to render the quantities in brackets domain, yields: 


nondimensional. Herein, b n is the half-chord of the 
airfoil and U n is the component of the flow speed, both 
normal to the elastic axis. 

In a simplified context, such a mixed form of the lift 
and moment was used in Simiu & Scanlan 9 , and 
Scanlan 8 . The unsteady aerodynamic derivatives for 
swept wings are obtainable from the previous equations 
of lift and aerodynamic moment, by assuming harmonic 




C La . ( 1 


a n \kj h - ib n -TjT 


tan A 


1 ^ - 1 ^ La„ . 

" 2 C L«„ a nfh +T — — a „ b n tan A 

“ * n 

(S r +\)^rik n + 5 r b n ?-J± tanA 
<?V dy 1 


time dependence of displacements quantities. In such a M = - 
way, the frequency domain counterpart of Eqs. (32)- 
(33), expressed in compact form, becomes: 


\ + Q n 


h L h {y,k n ,x)= pUlkX 
s M a {y,k n ,x)= pUlkX 


T~L\ +a,)L 2 

\ b - 






— M ] + a 0 M 2 
b„ 


(34) 

(35) 


tfcc 


if a "I ~~ a n \k n f a -ib n — tanA 


Jy 


, (40 a) 


j £ C La„ fa 


-(A-) 


J .r . L 2 

— f - r ~ a 

2 


2 k n 


if a + ^ ^La„ a n fa 2 k 2 ^ 


1 u $fa 
% 


Herein, the unsteady aerodynamic complex coefficients 
Li and Af, can be expressed as: 

L , = ///, + H a ; L 2 = /// 2 + // 3 ( 36 ) 


tan A 


(\ _ > 

1 C La n 1 


bj 

U n j 

l 2 k i i 


M x =iA x +A 4 \ M 2 = iA 2 + A 3 . (37) 


(«,+!)%-«*. +5 r *.^£- tan A 


(40 b) 
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Using in Eqs (39)-(40) equations (2), separating the real 
and the imaginary parts of the above expressions, the 
unsteady aerodynamic derivatives result as: 




H, = — 


tanA^' 
H k n dy 


£M + L {l+S f 

k„ 2 


G(k„) 


1 


+ b„ -^-tan A| 


d>’ 

«3 =-Q«, 

b n %a 


T-a* f( k n)fa +— 7^ fa + ~fa 




2 G(k„) 


F(k„), G{k „)( 1 _ ^ 




£ „ 

2 


[/a 2 $ & 


F ( k n\--<*n + T^r 


-| ( an 

kl dy 


tan^ 


2 * 2 <* 2 


«4=C t „ 


I A+ * 


c(*J 


/?„ tan A 


7 * 


/=■(*„) 


ff h 

dy 


+ 5, —b„ — ^-tan A 

2 <?y 2 


Y\ 


(41) 


A, — 


Q«, ( 


(i +a " 


tan A 


A/ 


- + 


tl) " 

*„ 2 rtW 


A? = 


C La. ^ 


2 k n 


° n 2 


1 

^ a n 

2 n 


2G(k„) 


fa 


-[al -^F{k n )f a +b n ^-xan\ 

A 3 = ^TT (0 + 2a„ M*. )/„ + ( 4 « 2 - l)c(*„ >„/„ 

4k; 

+ [ 2a " 

-(l-5 r 2 fl J)-<5A 2 f! + 2 a 0^tan 2 A) ’ 


Aa — — C r 


-<*nfh + - + *„ 


G(k n ) 


dy 1 

fh -trtanA 


(i + “-] f ' (t - ) t +i 'r A ^“ nA 




(42) 


The coupling terms, due to the sweep effects, are 
separated in the overall expressions of aerodynamic 
derivatives. For A = 0 the expressions of aerodynamic 
coefficients corresponding to straight wings are 
obtained. The unsteady aerodynamic derivatives 
recorded above coincide with the ones obtained 
differently in Barmby et al. 1 , where the span wise rates 
of change of bending and twist, a and X, that are 
associated with the sweep effect, also intervene. When 
specialized for 8 r =0, there coincide with the ones by 
Bisplinghoff et al. 2 

Aeroelastic Response of an Airfoil Featuring 
Plunging and Pitching Coupled Motions to Sonic- 
Boom Pressure Pulses. 


An applications on the aeroelastic response of an airfoil 
in an incompressible flow featuring plunging and 
pitching coupled motions to sonic-boom and blast 
pressure pulses will be given in the next developments. 
The aeroelastic governing system of equations of an 
airfoil featuring plunging and twisting degrees of 
freedom to sonic-boom, blast pressure pulses, expressed 
in dimensionless form, can be cast as: 

'(r)+ *a«'(0+ (5T/V t '(r) 

+ (5J/V , ) 2 ^(r)-/,(r) = /,(r) 

be a / r a )S '(r )+ a'(r)+ (2 C« l v )* '(* ) 

+a/v 2 -m a (r)= 0 ’ 


In the above expressions this following dimensionless 
parameters have been used: 

l„=L h b/mUl ; p m = P m b/rnUi ; V=Ujbco a ; 

H = m/xpb 2 ■ = cj2nm h ; m ; 

co h ={K h /m) n x co a =(K a /lJ 12 ; =cj2l a (O a ; 

r a = i/a / mbl )' 2 ’Z a =S a / mb ; m a = M a b 2 j I a U l . 


The dimensionless sonic-boom overpressure signature 
of the N-wave shock pulse, can be described as follow: 




■8 l ,H{x-rT p )o m 


/ A 

l-L 


•(45) 


Herein, the Heaviside step function H(t) has been 
introduced in order to describe the typical pressure 
time-history for blast or sonic-boom loads; <5^is a tracer 
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that should be taken as one when the sonic-boom is 
considered, and zero when the blast load is included; 
p m denotes the dimensionless peak reflected pressure 
in excess of the ambient one; T p denotes the positive 

phase duration of the pulse measured from the time of 
impact of the structure; r denotes the shock pulse 
length factor. For r = 1 the N-shaped pulse degenerates 
into a triangular pulse which corresponds to an 
explosive pulse (Fig. 4. a), and for r = 2 a symmetric 
N-shaped pulse is obtained. A depiction of l b jfp m vs. 
time is displayed in Fig. 4.b. 

The Eqs. (43) and (44) can be converted in the Laplace 
transformed space and solved for their unknowns, 

{=<£ (£)) and (a)); inverted back in time 

domain one obtain the plunging and pitching time- 
histories and the load factor time-history due to 

the sonic-boom pressure pulse, and 

{a(.y)}, respectively. 

When the dynamic response of the actively controlled 
lifting surface is analyzed, also the feedback control 
forces and moments, that are time dependent, have to 
be included in the Eqs. (43) and (44). This will be 
considered in the future work. 

RESULTS AND DISCUSSION 

Herein, an unified way enabling one to obtain the 
unsteady lift and aerodynamic moment in the time and 
frequency domains for swept aircraft wing was 
developed. This was done via the use of the indicial 
function approach. The time domain representation is 
essential towards determination of the dynamic 
aeroelastic response to time dependent external loads, 
and in the case of the application of a feedback control 
methodology, of the dynamic aeroelastic response to 
both external time dependent loads and control forces. 
The frequency domain representation is essential 
towards determination of the flutter instability. 

The unsteady aerodynamic derivatives for different 
values of a n , A and C La have been plotted in Figs. 5 

and 6, as a function of k n . For swept wings, the local 
lift-curve slope C La for sections normal to the elastic 

axis are obtained from the aerodynamics of swept 
wings (see Yates 10 ). 

The maximum influence of the corrective term 
(identified by the tracer <5 r ) is present the first plunging 
coefficients H { where the aerodynamic coefficient 
changes also its sign. Usually, for all coefficients, the 
effect of these terms becomes higher for high sweep 
angles. 


For swept wings, the local lift-curve slope C Lan for 
sections normal to the elastic axis are expressed as: 

C La n =Q«/ COsA - (49) 

The variation of the unsteady aerodynamic derivatives 
Hj and A i as a function of k n are depicted. In these 
developments, all the terms, including the aerodynamic 
ones associated with h and a , usually neglected, have 
been retained. As a result, the coefficients H 5 , H 6 and 
\ are also included. Whereas the aerodynamic 

coefficients //, and A 2 are the principal uncoupled 
aerodynamic damping coefficients in plunging and 
torsion, respectively, H 2 and A, are the coupled 
damping coefficients. For straight wings, these terms 
remains negative for all values of 2n/k n and for 

different values of a n , but for swept wings only A 2 
continues to remain negative. The elastic axis position 
is not involved in the expression of the unsteady 

aerodynamic coefficients //, and // 4 , a fact, which 
clearly appears from the equations. The variation of the 
aerodynamic derivatives as a function of k n is, in 
general, a smooth one. Among these coefficients, only 
H 2 features, with the variation of the reduced 
frequency k n , a change of sign. As concerns, the 

depiction of H t and A , versus 2 n/k n (= Ujnb n ), this 
representation enables one to get an idea of the 
variation of the respective quantity with that of the 
normal freestream speed U n . The differences in the 

unsteady aerodynamic coefficients induced by the 
discard or inclusion of the terms generated from the 
downwash velocity (i.e. that underscored by a solid line 
in Eq. (12)) are properly indicated in Figs. 7. As shown, 
the corrective term does not modify the trend of 

coefficients A x and A 4 . It should be mentioned that the 
expressions of the lift and aerodynamic moment in the 
frequency domain obtained by Barmby, Cunningham 
and Garrick 1 coincide with the ones obtained here via 
indicial function approach. 

The graphs depicting the aeroelastic response time- 
history to blast pulses (i.e. explosive and sonic-boom 
blasts) are displayed (Figs. 8 - 10). The graphs supply 
the dimensionless plunging displacement (^ = h/b), 
and the load factor (N = h*/g , where g is the 
acceleration of gravity). The predictions of £ and N 
based on pure plunging and coupled plunging-twist 
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models are depicted on the same plots, and the 
closeness of the two predictions becomes apparent from 
the graphs. Herein, the dotted and solid curves 
correspond to pure plunging and coupled plunging- 
pitch models, respectively. As a result, the coupling 
helps to reduce the amplitude of the aeroelastic 
response. The pitching displacement has its maximum 
for r= 0 The same trend is valid also for the load 
factor N. It should be indicated that the response to 
sonic-boom pressure pulse involves two different 
regimes; one for which 0 < r < 30” that corresponds to 
the forced motion, and the other one to T > 30” 
belonging to the free motion. The jump in the time- 
history of N is due to the discontinuity in the load 
occurring at r= 30”. This jump doesn’t appear for 
explosive pressure pulses, where r = 1. The increase of 
the mass ratio results in the increase of the plunging 
displacement amplitude and the decrease of the pitching 
amplitude. At the same time, for higher mass ratios, the 
differences in the plunging predictions based upon 
1DOF and 2 DOF disappear (Figs. 8 and 9). Moreover, 
for higher mass ratios, the motion damps out at larger 
times. Figs. 10 highlights the effect of the structural 
damping coefficient in plunging and pitching. 

Using the idea developed by Yates , a modified strip 
theory can be accommodated as to address the problem 
of the aeroelastic response (for open/closed loop 
aeroelastic systems), by capturing also the 3-D effects. 
Alternatively, an exact solution methodology enabling 
one to determine both the flutter instability and the 
aeroelastic dynamic response based on a double 
Laplace transform, in time and space, can be used. Such 
a method was devoted to the solution of aeroelastic 
eigenvalue problems in Karpouzian & Librescu 5,6 , 
Librescu & Thangjitham 7 . On this basis, the obtained 
results can be extended as to approach the subcritical 
aeroelastic response and flutter, respectively, of 3-D 
advanced lifting surfaces, in various flight speed 
regimes. 

CONCLUSIONS 

A unified treatment of the aeroelasticity of 2-D lifting 
surfaces in time and frequency domains has been 
presented and the usefulness in this context of the 
aerodynamic indicial functions concept was 
emphasized. 

Applications assessing the versatility of this approach 
enabling one to treat both subcritical aeroelastic 
responses and flutter instability were presented, 
and prospects for extending this treatment to 3-D 
aeroelastic problems are contemplated in forthcoming 
developments. 
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N(=h"/g) 



Fig. 7 Results based on the discard of the corrective 

effect (8 r = 0) , and including that effect (8 r = 1) 

— , vs. 2n/k n . 


pressure pulse 

H= 10, 1DOF; — = 10, 2 DOF; ^ = 50, IDOF; n = 50, 2DOF, 

{Pm = 1 ; ip = 1 5”; £h = £a = 0; V =1; r= 1; CO = 0.5; X a = 0; r u = 0.5; 
a = -0.2). 



Fig. 8 Influence of mass ratio on the time-history 
response of (a) the nondimensional plunging 
displacement and (b) of the load factor, to sonic-hoom 
pressure pulse 

H = 10, IDOF, jj. = 10, 2DOF; — [i = 50, IDOF; \i * 50, 2DOF, 

<£> m = 1;tp= 15”;Ch = ^a = 0; V =1; r = 2; <D = 0.5; Xa = 0; r« = 0.5; 
a = -0.2). 

11 


Fig. 10 Influence of structural damping coefficient on 
the time-history response of (a) the nondimensional 
plunging displacement, (b) of the load factor and (c) 
of the nondimensional pitching displacement to sonic- 
boom pressure pulse 

C h = Ca = 0;-.- Ch = C. = 0.1;. ^ = ^=0.25 Ch = C« = 0.5; 

(Pm» 15”; ? h = Co = 0; V =1; r = 2; CO = 0.5; Xa = 0; r a = 0.5; 
a = -0.2). 
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